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Search for Resonant Absorption of Solar Axions Emitted in Ml 
Transition in ^'^Fe Nuclei. 
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Abstract. A search for resonant absorbtion of 14.4 keV solar axions by ^^Fe target was performed. The 
Si(Li) detector placed inside the low-background setup was used to detect the 7-quanta appearing in the 
deexcitation of 14.4 keV nuclear level: j4-|-^'^Fe ^^Fe* ^^Fe-l-7. The new upper limit for the hadronic 
axion mass have been obtained: tua < 151 eV (90% C.L.) (5=0.5, 2=0.56). 
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1 Introduction 

A natural solution of the strong CP problem based on 
the global chiral symmetry U(l) was proposed by Peccei 
and Quinn [T]. Weinberg [2] and Wilczek [3] noted that 
spontaneous breaking of the PQ-symmetry at the energy 
/a leads to the existence of a new neutral spin-zero pseu- 
doscalar particle - axion. The axion mass (rriyi) and the 
strength of an axion coupling to an electron (gAe), a pho- 
ton (gA-y) and nucleons (gAN) are proportional to the in- 
verse of /a- The original PQWW-model with /a fixed at 
the electroweak scale /a = {2Gf)~^^^, was excluded after 
intensive experimental searches which were performed us- 
ing radioactive sources, reactors and accelerators (see [4] 
and refs therein). 

New axion models decoupled the scale of PQ symme- 
try breaking from the electroweak scale and the value of 
/a can be extended up to the Plank mass wlO^® GeV. 
Two classes of models for the "invisible" axion have been 
developed: the KSVZ or hadronic axion model [5j, 6J and 
the DFSZ or GUT axion models [T],^. 

The axion mass in both models is given in terms of 
neutral pion properties: 



rriA = — — ( 



Nl/2 



fA '{l + z + w){l + z)- 



(1) 



where /jr = 93 MeV is the pion decay constant, z = 
rriu/md ~ 0.56 and w = niu/ms ~ 0.029 are quark-mass 
ratios. For the given values of z and w the equation ^ 
is numerically presented as niA = 6.0 • 10®//a where uia 
and /a are in eV and GeV units, correspondingly. 

In contrast to the DFSZ axions, the KSVZ axions have 
no coupling to leptons and ordinary quarks at the tree 
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level, which results in the strong suppression of the in- 
teraction of the KSVZ axion with electrons through ra- 
diatively induced coupling Moreover, in some models 
axion photon coupling may differ from the original DFSZ 
or KSVZ gA'j couplings by a factor less than 10~^ [lOj . 

The axion coupling constants are constrained by vari- 
ous experiments, astrophysical and cosmological arguments. 
The astrophysical limits based on the axion interaction 
strength with photons and electrons in stars suggest that 
niA < 0.01 eV if one assumes the standard gA-y,Ae — "t-a 
relations (see ^JTV and refs therein). 

The axion-nucleon coupling is constrained by the up- 
per and lower limits based on the observed neutrino sig- 
nal of SN1987A [T2]. The supernova data leave open the so 
called "hadronic axion window" of ttia ''^ 10 eV if gAj cou- 
pling is sufficiently small. However, the cosmological limits 
on hot dark matter consisting of axions provide niA <1 
eV [13]. 

If axions do exist, the Sun would be an intense source 
of these particles. Several main mechanisms of solar axion 
production are considered. Axions can be efficiently pro- 
duced by Primakoff conversion of photons in the electric 
field of plasma. The resulting axion flux has an average 
energy of about 4 keV and can be detected by inverse 
Primakoff conversion in laboratory magnetic fields [14] - 
[17] or by the coherent conversion to photons in crystal 
detectors [l8]-[2T]. The experiments are sensitive to (7A7 
coupling. The obtained upper limits for the axion-photon 
coupling constant are gA-y < 10^^" 10^^ GeV"^, which 
still corresponds to the immense estimated axion flux at 
the level of 10"-hl0i3 cm-^s-^keV"!. 

The other source of axions are reactions of the solar 
cycle that produce solar energy [22] . Since axions are pseu- 
doscalar particles they can be emitted in nuclear magnetic 
transitions. Attempts to detect 478 keV monochromatic 



2 



A.V. Derbin, A.I. Egorov et al.: Search for the resonant absorption of solar axions. 



axions through the resonant absorption in '^hi nuclei tar- 
get have been performed in f23]-[26j. 

Since the temperature in the center of the Sun is — 1.3 
keV and some nuclei having low- lying nuclear levels can 
be excited thermally [21]. Monochromatic axions can be 
emitted in the magnetic nuclear transitions from the first 
thermally excited level to the ground state. The aim of 
the present work is to search for the 14.4 keV solar axions 
emitted by the Ml-transition in "''^Fe nuclei. The axions 
on the Earth can be observed in the inverse reaction of 
the resonant absorption via the registration of 7-rays (or 
conversion electrons) produced by the discharge of the ex- 
cited nuclear level [28]. The probability of emission and 
subsequent absorption of the axion in a magnetic tran- 
sition is determined only by the axion-nucleon coupling. 
The previous searches for 14.4 keV axion were performed 
in [2S1-1S2I. 



2 Emission and absorption of axions in 
nuclear transitions of magnetic type. 

As it has been found in [27], the most intense solar axion 
flux is connected with Ml-transition of ^^Fe nucleus. The 
energy of the first excited nuclear level 3/2 ~ is equal to 
14.413 keV, and the admixture of the E2 transition is 5^ = 
0.22%. The total axion flux <Pa depends on the level energy 
E^= 14.413 keV, temperature T, nuclear level lifetime 
= 1.34 /is, the abundance of the ^^Fe isotope on the Sun 
TV and the branching ratio of axions to photons emission 



by the expression [35 j l36 l l27] : 



uja/coj mm 



N 2exp{-E^/kT) 



(1 -I- 2exp{-E^/kT)) io-y 



(2) 



Owing to the Doppler broadening, the axion spectrum 
is a sum of Gaussian curves <Pa{Ea) with the dispersion 
cr(T) = E^{kT/My^'^, where T is the temperature at the 
point where the axion is emitted and M is the mass of 
the ^^Fe nucleus. The axion flux was calculated for the 
temperature dependence on the radius given by BS05(OP) 
Standard Solar Model [33] based on the corona high-Z 
abundances [51]. At the Earth, the differential axion flux 
at the maximum of the quasi-Gaussian distribution can 
be presented as: 



$a{Emi) = 4.15 X 10^^ ( — ) cm"2g-i]^gY-i_ 



(3) 



The obtained width of the axion line is equal to 175 = 2.2 
eV. This value significantly exceeds the energy of recoil- 
nucleus (1.8 /xeV), as well as the Doppler broadening of 
the line at temperature T = 300 K of the target nuclei 
(10 meV), and the self width of the level E = 4.65 • lO'^ 
eV. Thus, the percentage of axions satisfying the resonant- 
absorption condition amounts to the value ~ E / as- 

Within the framework of the long-wavelength approx- 
imation, the axion emission probability {uja/^-i) is given 



LOA_ 
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where Pj and pA are the photon and axion momenta re- 
spectively; a « 1/137, ^0 ~ 0.88, /is « 4.71 are the 
isoscalar and isovector nuclear magnetic momenta, /3 and 
r] are the parameters depending on the nuclear matrix el- 
ements. The values /? = —1.19 and 77 = 0.8 for the Ml 
transition in the ^^Fe nucleus were calculated in [27] . 

The interaction of the axion with nucleons is deter- 
mined by the coupling constant gAN, which consists of 
isoscalar g'^j^ and isovector parts. In the hadronic 
axion models, gAjsf and gj^jy constants can be represented 
in the form 0W\ : 



9% - -^[2S + i3E - D)^- 

GjA 1 



2w^ 



9\n = -^[{D + E) 



1- z 



(5) 



(6) 



2/a" ' ' '1 + z + w' 
Here, mjv = 939 MeV is the nuclcon mass, the constants 
D and E are expressed in terms of the isovector (Eas) 
and isoscalar (Eaq) pion-nucleon coupling constants. The 
exact values of D and E parameters determined from the 
semileptonic hyperon decays are equal to D = 0.808 ± 
0.006 and E = 0.462 ± 0.011 [37]. 

The parameter S characterizing the flavor singlet cou- 
pling still remains a poorly constrained one. Its value varies 
from S ~ 0.68 in the naive quark model down to 5 = 
—0.09 which is given on the basis of the EMC collabo- 
ration measurements |38) . The more stringent boundaries 
(0.37 < S < 0.53) and (0.15 < S < 0.5) were found in 
[39] and [40] , accordingly. As a result the value of the sum 
(SAAr/^ + siw) ™ O ™ay significantly decreases and, due 
to negativity of the parameter /3, actually vanishes. Tak- 
ing into account that the usually accepted value of u- and 
d-quark mass ratio z=0.56 can vary in 0.35-^0.6 range [4], 
the exact interpretation of experimental results is signifi- 
cantly restricted. We use S — 0.5 and 2=0.56 as reference 
when calculating the axion fiux for the KSVZ axion model. 

According to ([T]), the constants 3^4^ and gAN can be 
expressed in terms of the axion mass (S = 0.5) as 

(7) 
(8) 

The values of and gj^jy the DFSZ axion de- 
pend on the additional unknown parameter cos^ f3 which 
is defined by the ratio of the Higgs vacuum expectation 
values, but they have the same order of magnitude [9lll0j. 

Because the branching ratio (uJA/^-y) is a model depen- 
dent we can consider the parameter (ffl47v/5 + 5iAr)^ in the 
expression Q as a free unknown parameter characterizing 
the axion-nucleon coupling. 

The cross-section of the resonant absorption of the ax- 
ions is given by the expression similar to the one for the 
7-ray absorption and corrected by the loaI'^i ratio |41| : 



g% = -4.03 • 10-«(m^/l eV), 
g\^ = -2.75 ■ 10-"(m^/l eV). 



a{EA) — 2\fnaQ^ exp 



p2 



LOA_ 

^7. 



(9) 
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where (T07 is the maximum cross-section of the 7-ray res- 
onant absorption and F = l/r. The experimentally ob- 
tained value of (To-y for ^^Fe nucleus is equal to 2.56 x 10^^* 
cm^. 

In order to obtain the total cross-section value, we 
should integrate the expression ([9]) over the axion spec- 
trum given by ([3|). The integration of the narrow Gaus- 
sian distribution ^ over the wide axion spectrum yields 
a value close to <l>A{E]\fi). Using the dependence of <I>a © 
and (j{Ea) ^ on the {uiA/^-y) ratio, and, therefore, on the 
axion mass niA (|4l7l8p . one can numerically present the 
estimated rate of resonant absorption of axions by 5'^Fe 
nucleus {S=0.5, z=0.56): 



R = 1.56- 10~^(cjaM)^ 

5.16-m-Hg%l3 + gl^)HpA/p-,)' 
= 9.29 -10-34(^^)4(^^/^^)6 



(10) 

(11) 
(12) 



The amount of the detected 7-rays following the axion 
absorption is determined by the target mass, the time of 
measurement and the detector efficiency, while the obser- 
vation probability for the 14.4 keV peak depends on the 
background level of the experimental setup. 



3 Experimental setup. 

The planar Si(Li) detector with the diameter of the sen- 
sitive region 17 mm and 2.5 mm thick was used for the 
detection of 14.4 keV 7-rays. The detector was mounted 
inside the vacuum cryostat with 20 /im beryllium window. 
The 290 mg iron target enriched to 91% of ^^Fe isotope 
was placed directly on the beryllium window, the distance 
between the detector's surface and Fe target was «3 mm. 
The surface density of the target was 92 mg/cm^, while 
the attenuation of 14.4 keV 7-rays in iron corresponds to 
16 mg/cm^. 

The detector was surrounded by 12.5 cm copper and 
2.5 cm lead shields to eliminate the external 7-radioactivity. 
The background level at 14.4 keV was decreased in 110 
times in comparison with the unshielded detector. The 
detector was located above the ground surface, so in or- 
der to minimize the influence of cosmic radiation and fast 
neutrons, we used the active shielding assembled of five 
50 X 50 X 12 cm organic crystal scintillators. The neutron 
shielding and active muon veto cover the top and sides 
of the passive shield, except for the side where the detec- 
tor Dewar is located. The anticoincidence veto signal was 
obtained from the logical OR of all the photomultiplier 
tube discriminator outputs. The Si(Li) detector operated 
in an anticoincidence mode with plastic scintillators and 
the background events around 14.4 keV were reduced by 
a factor of 2.5. The rate of 50 /iS veto signals was 600 
counts/s, that leads to «3% dead time. The spectrum of 
the Si(Li) signals obtained in the coincidence with veto 
signals allows checking the probability of excitation 14.4 
keV level by the nuclear active component and cosmic ray 
muons. 

The spectrometric channel of the Si (Li) detector was 
organized in the following way. The first-stage field-effect 
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Fig. 1. Energy spectrum of the Si(Li) detector measured in 
the region 1-60 keV (main) and 7-20 keV (inset). 



transistor was mounted on a Teflon block, a few mm from 
the center contact of the silicon crystal, while the pream- 
plifier was placed beyond the passive shielding. The data 
acquisition system was based on standard Camac electron- 
ics. A signal produced by the preamplifier was transmitted 
to two separate amplifiers with different gain ratios, thus 
providing the possibility to collect spectra from both lower 
(0-60 keV) and higher energy (0-500 keV) regions. The 
amplifier outputs were converted using analog-to-digital 
converters, controlled by a PC through parallel interfaces. 
Taking into account the active shielding coincidence spec- 
tra for each amplifier, there were four 4096-channel spectra 
being recorded. 

The energy scale was defined using standard calibra- 
tion sources of ^^Fe, ^^Co and ^4iAm. The energy resolu- 
tion of the detector (FWHM) determined by the 14.4 keV 
7-line of ^^Co source turned out to be 280 eV. 

The detection efficiency of the Si(Li) detector was mea- 
sured with 14.4 keV 7-rays from a standardized ^^Co source. 
The self-absorption of 14.4 keV 7-rays by the iron tar- 
get was found via detailed M-C simulation. The over- 
all detection efficiency for 14.4 keV 7 is estimated to be 
(2.30±0.1)%. 



4 Results. 

The measurements were carried out during 84.5 days of 
livetime by ^2-hour runs. The obtained energy spectra 
are given in FiglT] Since the ^''Fe target consisted of a 
small amount of ^^Fe, the most intense peaks are con- 
nected with the characteristic X-rays of ^^Mn {Ka— 5.9 
keV, K/S—6A9 keV). One can clearly identify several peaks 
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Fig. 2. Fitting results of the energy spectrum inside the 12.4- 
16.6 keV region. The location of the expected axion peak is 
denoted by arrow. 



degrees of freedom at the 12.4 - 16.6 keV region amounted 
to 250. 

The fitting result is given in FigO The minimum 
value corresponds to the nonphysical value of the 14.4 keV 
peak area S = —25 events. The upper limit for the amount 
of events inside the peak was found via the conventional 
approach: the dependence of on the peak area S was 
calculated for various values of S while the rest parameters 
were left unrestrained. Then the appearance probability of 
the given x^(5') value was found and the obtained function 
P{x'^iS)) was normalized to unity for the S > region. 
Thus, the upper limit appeared to be equal to Sum — 24 
events. For the rate of axion absorption R given by pO)) 
the expected number of registered 14.4 keV 7-quanta is: 

S = e-rj-N-T-R<Sum- (14) 

Here, the number of ^^Fe nuclei N = 2.78 • 10^^, the mea- 
surement time t = 7.30 • lO^s, 7-ray registration efficiency 
e — 2.3 • 10^^ and internal conversion ratio rj — 0.105. 

The relation limits possible values of the axion- 
nucleon couplings constants and axion mass. In accor- 
dance with equations ()10|) and on condition that (jja/p-j)^ — 
1 provided for ttia < 2 keV one can obtain: 

I - 1.19 • g% + g\j^\ < 3.12 • 10^^ and (15) 



related to the ^^^Am radioactivity. The low-energy region 
contains the L-series of characteristic X-rays caused by 
Am -> 234Np* _^ ^ decay. The 13.9 keV peak consists 
of two lines with the energies 13.946 keV (13%, Lai) and 
13.761 (1.4%, Lq2), the 17.8 keV peak is more complex 
one formed by L^i_5 lines. The peaks with the energies 
of 10.84 keV and 13.2 keV (Ld and Lpi^^ of Bi, respec- 
tively) are present due to ^^^Pb— s-^^^Bi /3-decay in ^ssy 
series. Intense Am (59.5 keV) and ^lOfii (46.5 keV) 7- 
rays are observed at the higher energy region. All these 
peaks were used to determine the final energy scale and 
energy resolution ct(E) of the detector. 

FigO shows the detailed energy spectrum within the 
12.2 - 16.6 keV interval, where the axion peak was ex- 
pected. Apparently, there is no pronounced peak at the 
14.4 keV. In order to determine the upper limit for the 
number of events inside the expected peak we used the 
maximum likelihood approach. The likelihood function 
was determined as a sum of four Gaussians and the poly- 
nomial background. Three Gaussians represent the known 
characteristic X-rays, one gaussian stands for the expected 
14.4 keV axion peak and the third order polynomial is used 
for the smooth background: 



ruA < 151 eV at 90% c.l. 



(16) 
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Peak positions and energy resolution (a) were fixed, while 
the peak areas and background polynomial coefficients 
were independent free parameters. The total number of 



The limit ([TB)) is the strongest up-to-date result ob- 
tained with 14.4 keV solar axions. The previous limit {ruA < 
216 eV [31 ) is improved in 1.4 time. Because the expected 
intensity of 14.4 peak depends on m\ that in fact corre- 
sponds to an increase the sensitivity of the experiment 
in 4 times. The obtained limit on hadronic axion mass 
depends on the exact values of the parameters S and z 
(FiglH]). The uncertainty of the flavor-singlet axial- vector 
matrix element S changes the obtained constraints signif- 
icantly: ruA <217 eV (5=0.4) and tua <98 eV (5=0.68). 
Moreover, if the value of S is close to 0.17, the limit on the 
hadronic axion mass can not be derived from the present 
experiment. 

The value of u- and d quark-mass ratio 2=0.56 is gen- 
erally accepted for axion papers, but it could vary in the 
range 0.35-^0.6 [4]. For the fixed 5=0.5, the changes of z 
from 0.6 to 0.45 lead to the limits <(132-=-272) eV cor- 
respondingly. If z =0.34, the hadronic axion mass can not 
be restricted by the experiment. The reason is that the ra- 
tio {iJA/^-y) vanishes for some values of z and S. One can 
obtain from ([4]) and ([5]) that if S and z obey the relation 
5=1 — 1.5 • {z ± 0.01), the limit on the axion mass can 
not be obtained. Nevertheless, the limit on axion-nucleon 
coupling given by (|T5l) is still valid. 

As mentioned above, the main disadvantage of the ap- 
proach with axions emitted in 14.4 keV Ml-transition of 
^^Fe is that the nuclear-structure-dependent parameter (3 
has a negative value which, together with a poorly con- 
strained flavor singlet axial- vector matrix element 5, leads 
to large uncertainty of the (w^i/u;^) ratio. 
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Fig. 3. The hmit on the hadronic axion mass versus S for 
various values of z. The solid line corresponds to 2=0.56. 



5 Conclusion. 

A search for resonant absorption of the solar axion by 
^^Fe nuclei was performed using the planar Si(Li) detector 
installed inside the low-background setup. The intensity of 
the 14.4 keV peak measured for 84.5 days turned out to be 
< 0.28 events/day. This allowed us to set the new upper 
limit on the hadronic axion mass of tha < 151 eV (90% 
c.l.) with the generally accepted values S'=0.5 and z=0.56. 
The obtained limit strongly depends on the exact values 
of the parameters S and z. 
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